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Abstract 

This paper presents existence and uniqueness results for a class of parabolic systems with non linear diffu¬ 
sion and nonlocal interaction. These systems can be viewed as regular perturbations of Wasserstein gradient 
flows. Here we extend results known in the periodic case (0) to the whole space and on a smooth bounded 
domain. Existence is obtained using a semi-implicit Jordan-Kinderlehrer-Otto scheme and uniqueness follows 
from a displacement convexity argument. 


1 Introduction. 

In this paper we study existence and uniqueness of solutions for systems of the form 

f d t pi - div(/ 9 jV(V[p])) + a l div( / 9 1 Vy j '(/3i)) = 0 on R+x O, 

\ Pi{ 0, •) = Pi ,o on Q, 

where i £ [1, /] (7 £ N*), fl = R” or is a bounded set of R n and p := (pi,..., pi) is a collection of densities. Our 
motivation for this system comes from its appearance in modeling interacting species. 

In the case of V(V)[p]) = 0 or V, [p] does not depend of p, this system can be seen as a gradient flow in the 
product Wasserstein space i.e VU/ (pi) can be seen as the first variation of a functional T r defined on measures. 
This theory started with the work of Jordan, Kinderlehrer and Otto in EH where they discovered that the 
Fokker-Planck equation can be seen as the gradient flow of f Rn p log p + f Rn Vp. The method that they used 
to prove this result is often called JKO scheme. Now, it is well-known that the gradient flow method permits 
to prove the existence of solution under very weak assumptions on the initial condition for several evolution 
equations, such as the heat equation EH, the porous media equation [US, degenerate parabolic equations P, 
Keller-Segel equation [5]. The general theory of gradient flow has been very much developed and is detailed in 
the book of Ambrosio, Gigli and Savare, [2], which is the main reference in this domain. 

However, this method is very restrictive if we want to treat the case of systems with several interaction 
potentials. Indeed, Di Francesco and Fagioli show in the first part of [12] that we have to take the same (or pro¬ 
portional) interaction potentials, of the form V[p\ = W * p for all densities. They prove an existence/uniqueness 
result of m using gradient flow theory in a product Wasserstein space without diffusion (a* = 0) and with 
l = 2, Vi [pi, P 2 ] := W\,\ *pi + W \,2 * p 2 and V [pi , P 2 ] := V/ 2,2 * P 2 + W 2 P * pi where W \,2 and W 2 ,1 are equals or 
proportionals. Nevertheless in the second part of m, they introduce a new semi-implicit JKO scheme to treat 
the case where IT) .2 and W 24 are not proportional. In other words, they use the usual JKO scheme freezing 
the measure in Vi [p]. 

The purpose of this paper is to add a nonlinear diffusion in the system studied in 1 12 | . Unfortunately, this 
term requires strong convergence to pass to the limit. This can be obtained using an extension of Aubin-Lions 
lemma proved by Rossi and Savare in m and recalled in theorem 15.31 This theorem requires separately time- 
compactness and space compactness to obtain a strong convergence in L m ((0,T) x f 1). The time-compactness 
follows from classical estimate on the Wasserstein distance in the JKO scheme. The difficulty is to find the 
space-compactness. This problem has already been solved in 0 in the periodic case using the semi-implict 
scheme of [T2] . In this paper we extend this result on the whole space R™ or on a smooth bounded domain. On 
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the one hand in M™, we will use the same argument than in [ 8 j. We use the powerful flow interchange argument 
of Matthes, McCann and Savare m and also used in the work of Di Francesco and Matthes m- The differences 
with the periodic case are that functionals are not, a priori, bounded from below and we can not use Sobolev 
compactness embedding theorem. On the other hand in a bounded domain, the flow interchange argument is 
very restrictive because it forces us to work in a convex domain and to impose some boundary condition on 
Vi[p\. To avoid these assumptions, we establish a BV estimate to obtain compactness in space and then to find 
the strong convergence needed. 

The paper is composed of seven sections. In section [5] we start to recall some facts on the Wasserstein space 
and we state our main result, theorem 12.31 Sections [3l U and [5] are devoted to prove theorem 12.31 In section [3l 
we introduce a semi-implict JKO scheme, as in m, and resulting standard estimates. Then, in section 0] we 
recall the flow interchange theory developed in m and we find a stronger estimate on the solution’s gradient, 
which can be done by differentiating the energy along the heat flow. In section [5j we establish convergence 
results and we prove theorem 12.31 Section [G] deals with the case of a bounded domain. In the final section [71 we 
show uniqueness of using a displacement convexity argument. 


2 Wasserstein space and main result. 

Before stating the main theorem, we recall some facts on the Wasserstein distance. 


The Wasserstein distance. We introduce 

^(R™) := L Af(R n ;R + ) : f dp = 1 and M(p) := f \x\ 2 dp(x) < +oo 

l J R™ J K" 

and we note 7 ? f c (R n ) the subset of V 2 (R n ) of probability measures on 1 “ absolutely continuous with respect to 
the Lebesgue measure. The Wasserstein distance of order 2, W 2 {p, p), between p and p in V 2 (M. n ), is defined by 

W 2 (p,p):= inf (f \x-y\ 2 d^x^y)] 

7^n(p,/i) \J R n x R n J 

where II(p, p) is the set of probability measures on R" x R" whose first marginal is p and second marginal is p. 
It is well known that V 2 (R n ) equipped with W 2 defines a metric space (see for example CEiiiM). Moreover 
if p g V 2 c (M. n ) then W 2 {p, p) admits a unique optimal transport plan 7 t and this plan is induced by a transport 
map, i.e 7 t = (Id x T)#p, where T is the gradient of a convex function (see [B]). 

Now if p, p £ ^(K™) 1 ; we define the product distance by 

1 

W 2 {p , p) := ^2 W 2 (pi , pi), 

2=1 


or every equivalent metric as W 2 {p,p) := 1 W 2 (pi,pi). We can define also the 1-Wasserstein distance by 

Wi{p,p):= inf / \x-y\dj(x,y)) , 

7en(p,/i) \j K-xfn / 

and the Kantorovich duality formula (see uraim) gives 

Wi(p, p) = sup | J (fid(p-p) : ip £ ^(dlp - p\) n Lipi(R n )| , 

with Lipi is the set of 1-Lipschitz continuous functions. Then for all p, p £ V 2 Q& n ) and p £ Lip( R") we have 

[ Vd{p- p) ^CWi(p,p) ^CW 2 {p,p). (2.1) 

Jr™ 
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Main result. Let 1 £ N* and for all i £ [1,1], we define V) : V(W 1 ) 1 —> C 2 (R n ) continuous such that: 


• For all p = (pi, ...,pi)£ P(R n ) 1 , 

Vi[p\> 0, (2.2) 

• There exists C > 0 such that for all p £ P(R n ) 1 , x £ 1", 

l|V(Vi[p])|| L oo (H n) + \\D 2 (V[p})\\l~(m™) $ C, (2.3) 

i.e Vi[p\ and V(V)[p]) are Lipschitz functions and the Lipschitz constants do not depend on the measure. 


• There exists C > 0 such that for all v,cr £ ^(R™)*, 

IIV(VzM) - V(V5[<r])|| L . (R . ) < CW 2 (v,<t). 

Remark 2.1. The assumption (12.211 can be replaced by Vi[p] is bounded by below for all p. 


(2.4) 


Let m ^ 1, we define the class of functions by 

TL m := {x i-> xlog(x)} if m = 1, 

and, if to > 1, TL m is the class of strictly convex superlinear functions F : R + —> R which satisfy 

F(Q) = F'(0) = 0, F"(x) ^ Cx m ~ 2 and P{x) := xF'{x) — F{x) ^ C(x + x m ). (2.5) 

The two first assumptions imply that if to > 1 and F £ TL rn then F controls x m . 


Before giving a definition of solution of ED. we recall that the nonlinear diffusion term can be rewrite as 


div(pVF'(p)) = A P(p), 


where P(x) := xF'{x) — F(x) is the pressure associated to F. 

Definition 2.2. We say that (pi,...,p;) : [0,+oo[—» V£ c (M. n ) 1 is a weak solution of (11.11) if for all i £ [1,1], 
Pi £ C([0,T],P£ c (R")), Pi( Pi ) £ L 1 (]0,T[xR n ) for all T < oo and for all ip 1} ... ,ipi £ C c KJ ([0,+oo[xR’ 1 ), 


/»+00 r r- 

/ / [(d t pi - V<fi ■ \7(Vi[p})) pi + aiAipiPi(pi)] = - <pi(0,x)po,i(x). 

Jo J K" J K" 


With this definition of solution we have the following result 

Theorem 2.3. For all i £ [1,1], let Fi £ Ft mi , with rrii ^ 1, and V, satisfies El, El and (12.411 . Let 
ai,..., op positive constants. If po.i £ Vf c (R n ) satisfies 


Pi(po,i) + Vj(p 0 ,i|po) < Too, 


( 2 . 6 ) 


with 

■rt \ / F i(p( x )) dx if p <. C n , . | . f r i j 

™ := { Too otherwise, and V ^ := J R V ^]P d *- 

then there exist (pi,..., pi) : [0, +oo[— > Pf c (WL n ) 1 , continuous with respect to Wi, weak solution of (11.11) . 

Remark 2.4. In the following, to simplify the proof, we take on = 1. 
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3 Semi-implicit JKO scheme. 

In this section, we introduce the semi-implicit JKO scheme, as m, and we find the first estimates as in the 
usual JKO scheme. 

Let h > 0 be a time step, we construct l sequences with the following iterative discrete scheme: for all 
i £ [ 1 , Z], p ( j h = pifi and for all k ^ 1 , p k h minimizes 

^ApIpT 1 ) ■■= wliPiP^ 1 ) + 2 h {Mp) + MpIpT 1 )) » 

on p £ "P| C (R"), with p k h ~ l = {p\; h \ ..., pf^ 1 ). 

In the next proposition, we show that all these sequences are well defined. We start to prove that it is well 
defined for one step and after in remark f3.21 we extend the result for all k. 

Proposition 3.1. Let p 0 = (pip, ■. ■, pip) £ 'Pf c (]R n ) i ; there exists a unique solution p\ = {p\ h ,..., p] h ) £ 
Vf c {M. n ) 1 of the minimization problem above. 

Proof. First of all, we distinguish the case > 1 from rn l = 1. 

• If 77 ij > 1, then £i,h{p\p) ^ 0, for all p, pi,...,p; £ V% c (M. n ). Let p v be a minimizing sequence. As 
£i,h(pi,o\po) < +oo (according to (12.61) 1. (£i,h(pv\Po))v is bounded above. So there exists C > 0 such that 

0 < Ti(p v ) ^ C and W 2 (p v ,pi, o) ^ C. 

From the second inequality, it follows that the second moment of p„ is bounded. 

• Now if mi = 1, following m, we obtain 

£Mp) > \m{ P ) - C(1 + M(p)) a - \M{pl h ), (3.1) 

with some 0 < a < 1. And since x i —jX — C (l + x)“ is bounded below, we see that £ij t is bounded below. 

Let p„ be a minimizing sequence. Then we have (J-i(p „))„ bounded above. Indeed, as £i,h(Pi,o\Po) < +oo, 
{£i,h(Pu\Po))v is bounded above and from (12.21) we get, 

/ V i [p 0 ](x)p iy (x)dx > 0, 

Jr 

so {Fi{pv))v is bounded above. According to (13.11) . (M(p v )) v is bounded. Consequently (P,(p,f)) u is 
bounded because Ti(p) ^ —C( 1 + M(p))“. 


In both cases, using Dunford-Pettis’ theorem, we deduce that there exists p*^ £ Vf' (M n ) such that 

p v p\ h weakly in L 1 (M n ). 

It remains to prove that p) h is a solution for the minimization problem. But since T r and Wf (-, pip) are weakly 
lower semi-continuous in L 1 (R"), we have 

£ i,h{Pi,h\Po) < liminf £i,h(p v |p 0 )• 

v /'+oo 


To conclude the proof, we show that the minimizer is unique, 
and p £ P 2 C (® ra ) ^ Wi(p,plh) and the strict convexity of P L . 


This follows from the convexity of Vi(-|p 0 ) 

□ 


Remark 3.2. By induction, nronosition I H. 1\ is still true for all k ^ 1: 

the proof is similar when we take k — 1 instead of 0 and if we notice that for all i, 

?i{Pi,h) + Vi{pl h \pl) < ^i(Pi,o) + Vi(pi,o\p 0 ) + CW 2 ( Po , pi) < C. 
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The last inequality is obtained from the minimization scheme and from the assumptions m, m and HMD. 
By induction it becomes, for all k ^ 2, 

k -1 

Hpth 1 ) + Mpth 1 \p k h~ 1 ) < Mpi, o) + VifooiPo) + c'£ w M~ 1 ,pL) < c 1 - 

J =1 

27i*s inequality shows |p^ _1 ) < +oo and so we can bound (^(pu))^ in the previous proof. 


Thus we proved that sequences {p k h)k^s o are well defined for all i £ [1, ZJ. Then we define the interpolation 
Pi, h : R+ ->• V!f c (R n ) by, for all k £ N,’ 

Pi,h{t) = P%h if t£((k- 1 )h, kh]. (3.2) 

The following proposition shows that pi /, are solutions of a discrete approximation of the system JHD- 

Proposition 3.3. Let h > 0, for all T > 0, let N such that Nh = T and for all (</>i,..., fa) G C£°([0, T ) x R n ) z ; 
then 



Pi,h(t,x)d t fa(t,x) dxdt 


N-l 


k—0 

N-l 


-ft E / P l {p k i J f l {x))\(j)i{tk,x)dx 

fc=o , ' R " 

N-l 


E 

fc =0 


/ cc) dec, 

•/nr* 


with, for all (j> £ C£°([0,T) x R n ), 

\K[(j)\(x,y)\ < |||-D 2 </>|U- ([ 0 iT)xR r,)|a;- 2 /| 2 , 


and 7 ^ h is the optimal transport plan in T(p k h , p^ 1 ). 

Proof. We split the proof in two steps. We first compute the first variation of £i,h(-\p k ) and then we integrate 
in time. In the following, i is fixed in [1, ZJ. 


First step: For all k ^ 0, if 7*7 is the optimal transport plan in r (p k hi Pit, 1 ) then 


[ ViWiPifi'ix) ~Pi, h (x)) 
J R n 


h [ P i {p k i \ 1 {x))Ap i {x) dx 

J K" 

h [ S7{Vi[p k h ]){x) ■S7p i {x)p } l+ 1 {x)dx 
J R n 

[ P[ ( Pi\{x,y)d'y k h {x,y), 

JR 


R n x R n 


for all ifi G C£°(R n ) and for all i G [1, /]. 


~>oo /ttt) n irpn 


To obtain this equality, we compute the first variation of £i t h{-\p\). Let £ C° 
let T r defined by 

dpi T = & o T t , T 0 = Id. 

After we perturb p ^ 1 by p T = (T r )up^ 1 . According to the definition of p i ~l 1 , we get 


) and r > 0 and 


l^iAPr^-SiMp^U)) >o. 


(3.3) 


5 


By standard computations (see for instance |11| . |Tj) we have 


limsup -(W%(p T ,pl h )-Wiip^^lh)) ^ [ {y-x)-ii{y)d^ h {x,y), (3.4) 

T\0 T ' J^n x R n 


with 7 ^ is the optimal transport plan in W 2 (p^ h ,p i X ), 


limsup -(Fi{pr) ~ FiiPih 1 )) < - [ PiiPih 1 ^)) div (&0)) dx i 
t\ 0 T ’ J R" 


and 


If we combine (13.31) . (13.41) . (13.51) and (13.61) . we get 


(3.5) 


limsup - (Vi{p T \p%) - ViiPi^lPh)) < / V(Vi[Ph\)(x) ' Zi^PiVix) dx. (3.6) 

r\,0 T V / Jr™ 


[ (y-x)<i(y)d-/i }h {x,y) + h f \7iy i [p k h }){x)-£,i{x)p } lX 1 {x)dx-h ( P i {p’lX 1 (x)) div(&(x)) dx ^ 0. 

il" X E" JR" ’ J R" 


And if we replace by —this inequality becomes an equality. 

To conclude this first part, we choose ^ = \7ifi and we notice, using Taylor’s expansion, that 

<Pi{x) - tpi(y) = V<pi{y) -(x-y) +7Z[pi\{x,y), 


with 7 Z\ipi] satisfies 

\K-YPi\{x,y)\ < ^\\D 2 i Pi\\L°°([o,T)xR™)\x - y\ 2 ■ 

• Second step: let h > 0, for all T > 0, let N such that Nh = T ( tk '■= kh) and for all ((f>i,..., (f>i) G 
C“([0,T) x K") z , extend, for all i G [l,i], by </>i(0, •) on [—/i, 0), then 


/ Pi,h(t, x)d t (t>i(t, x) dxdt = V7 / p k h (x)d t 4>i(t 1 x)dxdt 
JM. n k —0 ^ ^k — i •'R 71, 

JV 

= E / Pi,h(, x )(Mtk,x) - </>i(t k -i,x))dx 


N—l 


E/ (f) i (t k ,x)(pl h (x)-p k i 1 (x))dx- Pi,o{x)<j)i(0, x) dx. 
_ q J R 71 J R 71 


Using the first part with (fi = (f>i(tk ,•), we get 


AT—1 


[ Pi,h(t, x)d t <f>i(t, x) dxdt = -/iV f P i (p k i X 1 (x))/±(j)i{t k ,x)dx 
J R 71 ./R 71 

AT—1 

7,._n R 71 


/c—0 
AT-1 


E 

/c—0 




/ p ifi (x)cl>i(Q,x)dx, 
J R" 


for all * G [1,ZJ. 
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□ 


The last proposition of this section gives usual estimates in gradient flow theory. 

Proposition 3.4. For all T < +oo and for all i £ [1,1], there exists a constant C < +oo such that for all 
k £ N and for all h with kh ^ T and let N = we have 


M(pl h ) < C, 

MpU < c, 

Ef=o 'wiip^p^^Ch. 


(3.7) 

(3.8) 

(3.9) 


Proof. The proof combines some techniques used in na et [12], In the following, i is fixed in [1,/]. As p i f l 1 is 
optimal and p\ h is admissible, we have 


^Ap:X\ph^ £ MM)- 


In other words, 


\wM, h ,p k c h l )+h (r^x 1 )+npU^ph)) < h (K(pth) + npiM)) ■ 


From (12.31) . we know that Vi[p\ is a C-Lipschitz function where C does not depend on the measure. Hence, 
because of (EH), we have 

Vi{pl h \p k h) - Viip^X 1 \p k h ) < CW 2 {p\%\pl h ). 

Using Young’s inequality, we obtain 

1 


WM) - Vi(Pi£ \Ph) < c*h + -wj( P % 1 ,Pi h ). 


It yields 


1 


- A wl(pt h , <r) < M-W,/.) - Wi, I 1 )) + 


(3.10) 


Summing over k, we can assert that 
^1 


/N -1 


Ei wZiplh’piX 1 ) < h E (^(' o a)-^(^,l 1 ))+ c ' 2r 


k =0 


, k—0 


< hiTM-T^ + ^T). 

But by assumption, JFi(pi, o) < +oc and —Ti{p) ^ C(1 + M{p)) a 1 with 0 < a < 1, then 
N 1 


E jWtipU’PiX 1 ) < h (■^(Pi.o) + C(! + M(p? h )) a + C 2 T ). 


(3.11) 


k =1 


Thus we are reduced to prove (13.71) . But 

M{p k i>h ) < 2IU 2 2 (p^,p i ,o) + 2M(p i) o) 


k -1 


^ 2 fc E + 2M (pi,o) 


m =0 


< 8 kh (Fiipifl) + C(1 + M(pl h )) a + C 2 T ) + 2 M(p it0 ) 

^ 8 T (Ftipifl) + C{ 1 + M(p* h ))“ + C 2 T) + 2M(pi, 0 ). 

As a < 1, we get (1X71) . The second line is obtained with the triangle inequality and Cauchy-Schwarz inequality 
while the third line is obtained because of (13.111) . So we have poved (13.71) and (13.91) . 

To have (13.81) . we just have to use (13.101) and to sum. This implies 

Fi(pi, h ) ^ FiiPifi)+C 2 T, 

which proves the proposition. 

□ 
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4 K-flows and gradient estimate. 

Estimates of proposition 13.41 permit to obtain weak convergence in L 1 (see proposition 15.ID . Unfortunately, it 
is not enough to pass to the limit in the nonlinear diffusion term Pi(p,;,h). In this section, we follow the general 

strategy developed in na and used in m and 0 to get an estimate on the gradient of ■ This estimate 
will be used in proposition 15.21 to have a strong convergence of p,;y,, in L mi (]0,T[x R n ). In the following, we are 
only interested by the case where rrii > 1 because if m; = 1 , Pi(pi,h) = Pi,h and the weak convergence is enough 
to pass to the limit in proposition 13.31 In the first part of this section, we recall the definition of k-Hows (or 
contractive gradient flow) and some results on the dissipation of J-, + Vi then, in the second part, we use these 
results with the heat flow to find an estimate on the gradient. 


4.1 k- flows. 

Definition 4.1. A semigroup © 5 , : R + xfj'fR") —»• p 2 C (K") a n-flow for the functional T : P 2 C (K") — > 
RU{+oo} with respect to W 2 if, for all p £ ^die curve s 1 —> ©(^[p] is absolutely continuous on R + and 

satisfies the evolution variational inequality (EVI) 

| CT=S w|(6^[p],p) + ^wi(GUp],p) < m (4.1) 

for all s > 0 and for all p £ P 2 C (R n ) such that H/(p) < + 00 , where 


d + ^ ,. f(t + s)-f(t ) 

-77 fit) ■= hm sup-. 

dz s—>-o+ ^ 


In [2], the authors showed that the fact a functional admits a K-flow is equivalent to A-displacement convexity 
(see section [7] for definition). 

The next two lemmas give results on the variations of p\ h along specific ft-flows and are extracted from m. 
The goal is to use them with the heat flow. 


Lemma 4.2. Let U/ : Pf c (R n ) — > RU{+oo} l.s.c on Pf c (R n ) which possesses a n-flot ©^. Define the dissipa¬ 
tion along ©^ by 

DiMn) := limsup - (P)(p) - J)(©^[p]) + Vj(p|/i) - Vi(6^[p]|/u)) 

«V> S 


for all p £ P 2 ac (R n ) and p. £ V% c (R n ) 1 . 

If fL~h Pi h are t wo consecutive steps of the semi-implicit JKO scheme, then 




k-l\ 


i , h ) - *(pth) > 


k -1 


) + $w2(pt h ,p*~ h 1 ). 


Proof. Since the result is trivial if T(pJ : ?l 1 ) = + 00 , we assume A>(p ih 1 ) < + 00 . 
Thus we can use the EVI inequality EH with p := p k ih and p := p i h x . We obtain 




1 d + 


J i,h ) ~ y 77 l°-=s ^2 (®l '\Pi,h\iPi,h ) + 2^ 2 \Pi,h\'Pi,h )• 

By lower semi-continuity of T, we have 

vKp*- 1 ) - Af(Pi, h ) > ^{P k ~h) - limiDf *{&%[pl h ]) 




limsup (^(pf- 1 ) - 4-(©i(p^))) 


^ lim sup 


1 d + 


*=. w^&UpiApiM 1 )) + 7^2 (pt h ,p*z 1 )- 


\,o v 2 do 


i,h 


k- 1 \ 


i,h 


(4.2) 
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The last line is obtained thanks to the WVcontinuity of s H 6 in s = 0. Moreover, the absolute continuity 
ofse> ©| ,\_Pih\ implies 


iim sup | ff=s W^GUpUlpth 1 )) > limsup -(w^&UpihiptH 1 ) Wiip^p^ 1 )) . 

s \0 V 2 da J a \0 5 V ' 


But since p\ h minimizes £i : h(-\Ph ), we get, for all s ^ 0 


wi(e\ 


[plhlpth) - w Z(plh,pth) > 2h PMh) - ^(6Uph J)) + 2 h {Vi{pl h \p k h - X ) - Vi{.<S%[pl h ]\p k h - 1 )) ■ 


This concludes the proof. 


□ 


Corollary 4.3. Under the same hypotheses as in lemmalet ©^ a n-flow such that, for all k G N, the curve 
s ©t,Lpf/J Has in L mi ( R n ), is differentiable for s > 0 and is continuous at s = 0. 

Let : Pf c (R n ) i—>■] — oo, +oo] be a functional such that 

iimmf (M6Upih])+n6Upth]\P k h - 1 ))) > ^ApUpV)- (4-3) 

Then, for all k G N, 


n^ 1 )-npi,h)>f^ i Mpth\p k h 


k -1 


) + 2 wM, h ,P k - h 1 )- 


(4.4) 


Proof. It is sufficient to show that Di^(-\p k *) is bounded below by Ai^(-\p k x ). The proof is as in corollary 
4.3 of P3J. The hypothese of L mi -regularity on ©^ imply that s K > ^(©'ij, [p k },]) is differentiable for s > 0 and 
continuous at s = 0. We have the same regularity for s Vi(©^[pf h ]|p^ _1 ). By the fundamental theorem of 
calculus, 


ViApihlpt 1 ) = hmsupi(J- i (p^)-J- i (©^[^]) + V i (p^|p^- 1 )-V i (6^[^]|^- 1 )) 

8\0 S 




v f 1 j 

( d 

lim sup / 

K du |„= 

0 

0 

/ 

05 

r 1 ,. r ( 

' d 

/ lim inf 


0 

/ 

cc 

0 

, dcr\„ = 


{Fi{&l[plh\) + Vi(.&l[Pi,h]\Ph~ 1 )) dz 


The last line is obtained by Fatou’s lemma and assumption (14.31) . To conclude we apply lemma 14.21 


□ 


4.2 Gradient estimate. 


Proposition 4.4. For all i G [1, i] such that irii > 1, there exists a constant C which depends only on pi$ such 
that 

II PT,h I L 2 ([0,T] ;«!(«")) ^ C{ 1 + T) 

for all T > 0. 

Before starting the proof of the proposition 14.41 we recall the definition of the Entropy functional , 

E{p) = [ p log p, for all p G P ac (R"). 

J R» 

We know that this functional possesses a K-flow, with n = 0 which is given by the heat semigroup (see for 
instance GO, S3] or [l5]). In other words, for a given 770 G Vf c (M. n ), the curve s 1 —> rj(s) := ©^[ 770 ] solves 

f d s rj = Arj on R + x R n , 

\ 77 ( 0 ) = 770 on R m , 

in the classical sense. 77 (s) is a positive density for all s > 0 and is continuously differentiable as a map from 
R + to C°° D T 1 (R"). Moreover, if 770 G L m (W l ), then rj(s) converges to 770 in L m (]R") when s \ 0. 
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Then we define 


Proof. Based on the facts set out above, &e satisfies the hypotheses of the corollary 14.31 We just have to define 
a suitable lower bound &j, E to use it. The spatial regularity of rj(s) for all s > 0 allows the following calculations. 
Thus for all p € V% c (W n ) 1 , we have 

d s (J r i (& s E [r]o])+ Vi(& s E [rio]\tx)) = [ d s Fi(r])dx+ [ Vi[p)d s ri(s,x) dx 

J R 71 JR 71 

= / El(rj(s,x))Ari(s,x) dx + / Vi[p]Ari(s,x) dx 

J R 71 J R 71 

= - [ E"(ri(s,x))\Wri(s,x)\ 2 dx + f A(Vi[n])ri(s,x)dx. 

J R 71 J R 71 

According to (12.51) . F"(x) > Cx mi ~ 2 thus 

9 S (Ji(6|,[77o]) + Vi{& s E [r]o]\n)) ^ -C f r?(s, x) mi ~ 2 \Vr](s, x)\ 2 dx + f A(Vi[n\)ri(s,x)dx 

J R" J K" 

< — C f \Vrj(s,x) m C 2 \ 2 dx + [ A(Vi[n})r](s,x) dx. 

J R" J K" 

&i iE (p\li) := C [ |V(p(s,x) mi/2 )| 2 dx - [ A(V\n\)p(s,x)dx. 

JM" J R" 

We shall now establish that J \^ E satisfies (14.311 . First of all, we notice that 

iimjnf (-^| _ (M&E[pih])+Vi( 6 E[pth}\P^ t )^ > (-^| _ WkiPi,^ ( 4 -5) 

+ 'rff - ( 4 . 6 ) 

Thanks to the proof of lemma 4.4 and with lemma A.l of [13], we obtain 

liminff-A (H6° E [pt h ))) [ \V(pt h (xr i/2 )\ 2 dx. (4.7) 

S A0 V d(7 U=» J J K" 

Moreover, as & S E is continuous in i 1 (K Tl ) at s = 0 and according to (12.311 . 

iiminf j^mP k h - l })pU x )dx. (4.8) 

The combination of (14.51) . (14.71) and (14.811 gives CO} for We apply corollarv l4.3l and we get 

E(P ^) ~ E(pZ h ) > hSUMPi^pt 1 )- (4-9) 

but since A(V)[p]) G L 00 (M rl ) uniformly on p (12.31) . 

Ch [ \V(pl h (xr i/2 )\ 2 dx < E{p\tf) - E{pl h ) + Ch. 

JR n 

Now we sum on k from 1 to TV = 

N 

Ch ||V(p^(x)^ /2 )|| 2 2(R „) < E( Pi , o) - E(p» h ) + CT. (4.10) 

k= 1 

According to HP and m, there exists a constant C > 0 and 0 < a < 1 such that for all p £ Vf c (M. n ), 

-C(l + M{p)) a ^E{p)^CFi{p). 

Since for all k, h, M(p l f h ) is bounded, according to (13.71) and the fact that J r i (p ii0 ) < +oo by (12.61) . we have 

N 

l|V(p.^(x) mi / 2 )|| 2 2(R „ ) < C(1 + T). 


k =1 


To conclude the proof, we use (12.51) and (13.811 . 


□ 
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5 


Passage to the limit. 

5.1 Weak and strong convergences. 

The first convergence result is obtained using the estimates on the distance ( 13 . 91 ) and on the energy T, ( 13 . 81 ) . 

Proposition 5.1. Every sequences (hk)ken of time steps which tends to 0 contains a subsequence, non-relabelled, 
such that Pi,h k converges, uniformly on compact time intervals, in W 2 to a ^-Holder function pi : [0,+oo[— > 
Vf c (R n ). 

Proof The estimation on the sum of distances gives us for all t,s < 0, 

W 2 {pi,h(t, ■), Pi,h(s, •)) < C(\t - s| + h) 1/2 , 


with C independ of h. 

According to the proposition 3.3.1 of [5] and using a diagonal argument, at least for a subsequence, for all 
i, Pi.h k converges uniformly on compact time intervals in W 2 to a ^-Holder function pi : [0,+oo[—>■ 

To conclude we show that for all t d 0, p(t, •) £ Pf c (R n ). But as Fi is superlinear, Dunford-Pettis’ theorem 
completes the proof. 

□ 


With the previous proposition, we can pass to the limit in the case mj = 1 because Pi(pi^h) = Pi,h and in the 
term V(V)[pi,fc,]) thanks to the hypothesis (12.41) . Unfortunately, it is not enough to pass to the limit in Pi(pi^h) 
when m-i > 1. In the next proposition, we use proposition 14.41 to get a stronger convergence. 

Proposition 5.2. For all i £ [Id] such that mi > 1, pi } h converges to pi in L mi (]0,r[xR“) and Pi(pi t h) 
converges to Pi{pi) in L 1 (]0,T[x R"), for all T > 0. 

The proof of this proposition is obtained by using an extention of Aubin-Lions lemma given by Rossi and 
Savare in |T7] (theorem 2) and recalled in [13| (theorem 4.9). 

Theorem 5.3 (th. 2 in |17|). On a Banach space X, let be given 

• a normal coercive integrand Q : X —> R + , i.e, G is l.s.c and its sublevels are relatively compact in X, 

• a pseudo-distance g : X x X —> [0, + 00 ], i.e, g is l.s.c and [g(p , p) = 0, p, p £ X with G(p), G{p) < 00 ] => 
p = p. 

Let U be a set of measurable functions u : ]0, T[— > X with a fixed T > 0. Under the hypotheses that 

r-T-h 

G(u(t))dt<+ 00 and lim sup / g(u{t + h),u(t)) dt = 0, (5.1) 

h \t°u£U Jo 


sup 

uet/ 


f 

Jo 


U contains a subsequence (u n )neN which converges in measure with respect to t €=]0, T[ to a limit it* : ]0, T[—>• X. 
To apply this theorem, we define on X := L mi (R”), as in 03, 9 by 


and Qi by 


f W 2 (p,p) if p,/i G P 2 (R”), 
+00 otherwise, 


Gi(p) ■■= 


||p mi / 2 ||i? l (R") + M(p) 
+00 


if p £ Vf c (R n ) and p mi / 2 S H\ R”), 
otherwise. 


Now, we show that G% satisfies theorem 15.31 conditions. 


Lemma 5.4. For all i £ [Id] such that mi > 1, Gi is l.s.c and its sublevels are relatively compact in L mi (J8L n ). 

Proof. The l.s.c of Gi on L mi (M. n ) follows from lemma A.l in p3J. To complete the proof we have to show that 
sublevels A c := {p £ L mi (M. n ) \ Gi(p) ^ c} of Gi are relatively compact in L mi (R n )- 

To do this, we prove that B c := {rj = p m G 2 | p £ A c } is relatively compact in L 2 (R”) and since the map 
j : L 2 ( R") —> L mi (R")> with j{rf) = g 2 / mi , is continuous, A c = j(B c ) will be relatively compact in (!(")■ 
We want to apply the Frechet-Kolmogorov theorem to show that B c is relatively compact in L 2 (R n ). 
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B c is bounded in L 2 (R n ): Since ry 2 = p mi with Gi (p) ^ c, it is straightforward to see 


/ 

J R r 


ry 2 < c. 


• B c is tight under translations: for every p £ B c and h G M” we have that 

f \p(x + h) — 77 ( 33 )| 2 dx < \h\ 2 f f \X7p(x + zh)\dz dx < \h\ 2 f \Xp(x)\ 2 dx ^ c\h\ 2 , 
JR" J R" JO JR" 

thus the left hand side converges to 0 uniformly on B c as \h\ \ 0. 

• Elements of B c are unifomly decaying at infinity: For all p £ B c and R > 0, we have 


[ dx < TTm [ |a:| 1/ "ry 1/ " m *»y 2 - 1/ " w *‘ dx. 
J\x\>R ft ' JR n 


' |x|>i? 

If we use Holder inequality with p = 2n and q = 2 n-i » we 

1/2 n 


J\x\>R U 1 \J R 


\x\ 2 p 2/r> 


j 

JR r 


2(2mi — l/n)/rrii(2—l/n ) 


As p 2 / mi = p with Qi(p) ^ c, we have 


[ \x\ 2 p 2/m ' < c. 
J R" 


To bound the other term we use the Gagliardo-Nirenberg inequality: for 1 ^ q, r ^ +oo, we have 


IMlLp^c'||v«||2 r |H| 


a in.nl—a 

5 


for all 0 < a < 1 and for p given by 


= a\ - 


1 1 


+ (l-a)-. 


We choose p = ^pf-ijn) ’ 9 = r = 2 and ck = 2 ( 2 ml-i/n) ( s i nce nii > 1 we have 0 < a < 1) then we 


obtain: 


f ^(2m i -l/n)/m i (2-l/n) ^ / f (V^V^ 2 ( f rf 
JR" VJr" / VJr" 

but since ry = p mi / 2 with Gi(p) ^ c, the second term is bounded then 


(l-a)p/2 


' |rc|>i? 


p 2 dx < 


C 

R l/n 


as R goes to +oo. 

We conclude thanks to Frechet-Kolmogorov theorem. □ 

Proof of the vrovosition \5.S[ We want to apply theorem [53] with X := L mi (R ra ), Q := Gi, g and U := {pi,h k \ k £ 
N}. According to lemma 331 Gi satisfies the hypotheses of the theorem. It’s obvious that it is the same for g. 
Thus we only have to check conditions for U. The first condition is satisfied because of (E7| and dm and the 
second is satisfied because of (13.91) (the proof is done in m proposition 4.8, for example). 

According to theorem 15.31 and using a diagonal argument, there exists a subsequence, not-relabeled, such 
that for all i with m; > 1, there exists pi :]0,T[— > L mi (M. n ) such that Pi t h k converges in measure with respect 
to t in L mi (M. n ) to pi. Moreover, as Pi,h k {t) converges in W 2 for all t £ [0,T] to pi(t) (proposition EH then 
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pi = Pi. Now since convergence in measure implies a.e convergence up to a subsequence, we may also assume 
that Pi t h k {t) converges strongly in L mi (R n ) to pi(t) t- a.e. Now, thanks to (13.81) and (12.51) we have 



dx ^ CJ'i(pi' h (t, •)) < C, 


then Lebesgue’s dominated convergence theorem implies that Pi t h k converges strongly in L mi (]0,T[x R") to pi. 

To conclude the proof we have to show that Pi(pi t h) converges to p i(pi) in i 1 (]0,T[x R n ). First of all, up 
to a subsequence, we may assume that there exists g £ L mi (]0,T[x R n ) such that 


Thus according to (12.51) 


Pi,h k -> Pi (t,x)-a.e and p i>hk < g (t,x)- a.e. 


Pi(pi,h k ) p i(Pi) (t,x)-a.e and 0 < P{pi,h k ) K C{Pi,h k + g mi ) (t,x)~ a.e. 


So when we pass to the limit we have (t, x)-a.e 


0 < P( P i) < C( Pl +g mi )£ L\] 0, T[x R"). 


Then C(p ithk + pi + 2 g mi ) - | p i{pi,h k ) ~ p i{Pi)\ > 0 and 


2CT + 2cff g{x) mi dxdt = ff liminf (C(p ijhk + Pi + 2g m *) - \Pi(pi, hk ) - Pi{pi)\) 

i/]0,T[xR" i/10,T(xR" 



JJ]0,T[x R" 



2CT + 2C 

If 

J ]0,T[x R’ 

< 

2CT + 2C 

ff 

J]0,T[x R’ 


/]0,T[x R 

[f 

U ]0,T[x I 


To do these computations, we used that \\pi,h k ||l 1 (]o ) t[x R") = l|Pi||L 1 (lo,T[x R") = p an d Fatou’s lemma. Since 
g £ L mi (]0,T[x R”), we obtain 


limsup // | p i(pi,h k ) ~ p i{Pi) I < 0, 

JJ]0,T[x R" 


which concludes the proof. 


□ 


5.2 Limit of the discrete system. 

In this section, we pass to the limit in the discrete system of propsosition 13.31 In the following, we consider 
<t>i £ C c °°([0,T) x R n ) and N = [£j. 

proof of theorem \2.S\ We will pass to the limit in all terms in proposition 13.31 


• Convergence of the remainder term: By definition of 1Z, we have 

[ p [Mtk, ■)}(x,y)d^ h (x,y) < ^||V 2 0, ||l~ ([ o,t ] xr-)W / 2 2 (p^, pff 1 ). 
J R n x R n z 

and according to the estimate (|3.9[h we get 


N—l 

E 


Tl[(j)i(t k , -)\{x,y)d^ h {x,y) 


1 X i 


k—0 

• Convergence of the linear term: 

r T r rT 


N—l 


o. 


k =0 


/ / Pi,h(t, x)d t (l)i(t, x)dxdt - / / pi(t,x)d t 4>i(t,x) dxdt 

Jo J R" Jo J R" 


< CT sup W 2 (pi,h(t, -),pi(t, •)) ->■ 0, 
te[o,T] 


when h \ 0 because of propostion 15.11 
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• Convergence of the diffusion term: 

N-l 


h y [ P i (p'l\ 1 (x)) ■ Acj)i(tk,x) dx - f f P i (p l {t,x))/S.(t) i {t,x)dxdt 
^,_q J R 71 JO J 

[ ( Pi(pi,h(t , x)) - Pi(p(t, x))) A 4>i(t, x) dxdt 
J K" 


^ CT\\D 3 (/)i\\L°°h + 


/ 0 J R' 


If to, = 1, the right hand side converges to 0 because of proposition 15.11 and otherwise it goes to 0 because 
of proposition 15.21 


• Convergence of the interaction term: 

N -1 




+ 


h'V] f S7 (Vi[p^\)(x) ■ V(/) i (tk,x)pi~l 1 (x) dx - f f V(Vi[p(t, •)])(£) • V0j(t, x)p i (t 1 x) dxdt 
_0 J R 71 J 0 J R 71 

-|V —1 « rtk+i r 

h J2 V(Vi[p%\)(x) ■ V(j) i {t k ,x)p^ 1 {x)dx -/ / V(T4[p(t,-)])(x) • V</> i (4,x)p' c i + 1 (x)dxdt 

i „—n JR 71 —n J tk JR 71 


k—0 

N ~ i ft k+ 1 


^2 [ f V(Vi[p(t, -)])(x) • (V(/>i(t fe ,x) - Vcj) i (t,x))p^ 1 (x)dxdt 
u—n Jtk JR 71 


k—0 

N-l ntk + l 


+ E 

fc=0 J 

^ J\ + J 2 + < 73 . 


V(Vi[p(i, •)])(a ; ) • V</> i (t,x)(p^ 1 (x) - pi(t,x))dxdt 


As converges weakly L 1 (]0,T[x K n ) to pi and (V(V*[p]) • V0j) £ L°°([0, T] x R ra ), then 

J 3 —^ 0 as h y 0 . 

For J 2 , we use the fact that Voh is a Lipschitz function and that V(U[p]) is bounded thanks to (12.31) . 
and then, 

J 2 ^ CT||L> 2 ^ i || i oo([ 0iT ] xR ^)/i —y 0. 

Using assumption (I2~l . we have 

r t k +1 

Jl ^ C f ||V0i||L°o([O,T]xR 71 ) ^ (Phi dt 

1 _n -/tfe 


fc=0 


^ 1 /*^fc+i 

< C||v&[| i « ([0i31x *» ) E / (^(pt p£ + 1 ) + w 2 (p£ +1 , p(t ,•)))dt 


fc=0 Jtk 


/ N 1 N 1 /»t fc+1 

< c'||v^|| i =o ([0 ,nxR») UE^Cp^p^^ + E / w 2 (p* + \p(t,-)) 

V fc =0 fe =0 ■'** 

< cuv^lUoo^TjxKn) (r]T w*(p k h ,p k h +1 ) +J*w 2 (p h (t,-),p(t,-))dt\ 


According to (13.91) . we obtain 


when h\ 0. Moreover, 


JV-l 


Tj2wi(ptp h h +1 )^lCTh^0 


k—0 


[ W 2 (p h (t, -),p{t, -))dt < T sup W 2 (p h (t, -),p{t, •)) -t 0, 
7o te[o,T] 

when h goes to 0, which proves that 


J 1 —y 0 as h —y 0. 

If we combine all these convergences, theorem 12.31 is proved. 


□ 
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6 The case of a bounded domain Q. 

In this section, we work on a smooth bounded domain LI of R n and only with one density but, as in the whole 
space, the result readily extends to systems. Our aim is to solve eh- We remark that f l is not taken convex so 
we can not use the flow interchange argument anymore because this argument uses the displacement convexity 
of the Entropy. Moreover since Ll is bounded, the solution has to satisfy some boundary conditions contrary to 
the periodic case [5] or in R ra . In our case, we study EH) with no flux boundary condition, which is the natural 
boundary condition for gradient flows, i.e we want to solve 

( d t p — div(pV(T[p])) — AP{p) = 0 on R + x fl, 

l (pV(V[p]) + VP{p)) • v = 0 on R + xdLl, (6.1) 

1 P(0>') = Po on R", 


where v is the outward unit normal to dL1. 

We say that p : [0, +oo[—>■ V ac {Ll) is a weak solutions of (16.11) . with F £ TL m , if p £ C([0, +oo[; V ac {Ll)) l~l 
L m {]0,T[xLl), P{p) £ ZdQO^xfJ), VP(p) e M n ([0,T] xLl) for all T < oo and if for all ip £ C“([0, +oo[x R"), 
we have 

[ [ [{dtp - V<p • V{V[p]))p - VP(p) ■ Vp\ = - [ p{0, x)po{x). 

J 0 J Q J 

Since test functions are in C(?°([0, +oo[x R ra ), we do not impose that they vanish on the boundary of f1, which 
give Neumann boundary condition. 


Theorem 6.1. Let F £ for m > 1 and let V satisfies : B . E3, (HU). If we assume that po £ V ac {{l) 
satisfies 


F{Po) + V(po|po) < +oo, 


with 


Hp) ■■= 


In F (p( x )) dx ifp<&£ [n, 
+oo otherwise, 


then EH) admits at least one weak solution. 


and V{p\p) := 



( 6 . 2 ) 


The proof of this theorem is different from the one on R n because we will not use the flow interchange 
argument of Matthes, McCann and Savare to find strong convergence since LI is not assumed convex. First, we 
will find an a.e equality using the first variation of energies in order to have a discrete equation, as in proposition 
13.111 Then, we will derive an new estimate on the gradient of some power of ph from this a.e equality. To conclude, 
we will use again the refined version of Aubin-Lions lemma of Rossi and Savare in |17| . 

On LI we can define, with the semi-implicit JKO scheme, the sequence {p^)k but this time we minimize 


p ^ Sniplpt 1 ) := ^W*{p, pi 1 ) + F{p) + Vfyl^" 1 ) 


on V{LV). The proof of existence and uniqueness of is the same as in proposition 13.II It is even easier because 
on a bounded domain T is bounded from below for all to ^ 1. We find also the same estimates than in the 
proposition 13.41 on the functional and the distance (see for example mm)- 


Now we will establish a discrete equation satisfied by the piecewise 
by, for all k £ N, 


mtcipuiaiiwn < 


1 \r ns 


Ph{t) = Ph, if t £ {{k - 1 )h,kh\. 


Proposition 6.2. For every k ^ 0, we have 

{y - T k {y))p k h +1 + hV{V[p k h })p k h +1 + hV{P{p k+1 )) = 0 a.e on LI, (6.3) 
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where Tk is the optimal transport map between and p\. Then ph satisfies 


N-l 


f f p h (t,x)d t <p(t,x) dxdt = h^2 [ ^( v [Ph))( x ) ■ ^^(t k ,x)p k h +1 (x)dxdt 

do J £1 _Q J 

N-l 

+h Tj j ^ P (Ph +1 ( x )) ' Vtp(t k ,x)dx (6.4) 

fe=o 
N-l r 

+ Tj / 'R[ i p(tk,-)\( x >y)d'Yk(x,y) 

k=0 Jnxn 

- / p o (x)ip(0,x)dx, 

J n 

with N = [^], for all </> £ C(?°([0, T) x R n ), 7 *, is the optimal transport plan in r(pj(, p j) +1 ) and 

\n[(p\{x,y)\ < i||£ ) 2 (/'||L-(RxR")|a; - y\ 2 - 

Proof. First, we prove the equality (16.31) . As in proposition 13.31 taking the first vartiation in the semi-implicit 
JKO scheme, we find for all £ £ C“(fl; R. n ), 

f (y- T k (y)) ■ ti(y)p k h +1 (y) dy + h [ v(v[p k ]) • f P k h +l - h [ p( P k+1 ) div(£) = o, (6.5) 

J £2 J J f2 

where Tk is the optimal transport map between p h +1 and p\. Now we claim that P(p^ +1 ) £ W /1,1 (fl). Indeed, 
since F controls x m and P is controlled by x m then (13.81) gives P{Ph +1 ) £ L 1 (f2). Moreover, (16.51) gives 


[ P(p k+1 ) div(0 
J o 




ly-Jfc^)! fc+i 

U Ph 


c 


IICIlL“(n) < 


W 2 (ptp k h +1 ) 


c 




This implies P{p k h +1 ) £ BV(Q.) and VP(p k+1 ) = VV[p k h \p k+1 + dd=2± p f f+ 1 in M"(!l). And, since S7V{p k h ]p k h +1 + 
Id f Tk p k+1 £ F 1 (fi), we have P(p k+1 ) £ W l d(Vf) and (16.31) . 

Now, we verify that ph statisfi.es (16.41) . We start to take the scalar product between (16.31) and V<p with 
p £ C£°([0, T ) x R n ), and we find, for all t £ [0,T), 

[ (y-Tk{y))-Vip{t,y)p k+1 {y)dy + h [ V(y[p k h \){y)-S7ip{t,y)p k h +1 (y)dy + h f V(P(p k+1 ))(y)-Vtp{t,y) dy = 0. 
J 17 J f2 J 

( 6 . 6 ) 

Moreover, if we extend ip by y>(0, •) on [—h, 0), then 


n 
1 0 


Ph{t, x)dtip(t, x) dxdt = / / Ph(x)dtip(t,x) dxdt 

k—0 ^ ^k-i J ^ 

N „ 

= / Phirfiviticx) - <p{tk-i,x))dx 

= T p(tk,x)(p k h( x ) - p k h +1 {x))dx~ / po(x)tp(0,x)dx. 
k=0 ^ ^ ^ K 71 

And using the second order Taylor-Lagrange formula, we find 

/ (ip(kh,x) - ip(kh,y))dj k {x,y) = Vtp(kh,y) ■ (x - y) dry k (x,y) + H[tp(t k ,-)](x,y)d'y k (x,y). 

JflxCl dflxfl «/f2xf2 

This concludes the proof if we sum on k and use m- 


□ 
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Remark 6.3. We remark that equality a is still true in R n . Indeed, the first part of the proof does not 
depend of the domain and we can use this argument on R™. This equality will he used in section 0 to obtain 
uniqueness result. 

In the next proposition, we propose an alternative argument to the flow interchange argument to get an 
estimate on the gradient of p^. Differences with the flow interchange argument are that we do not need to 
assume the space convexity and boundary condition on VP[p]. Moreover we do not obtain exactly the same 
estimate. Indeed, in proposition 14.41 Vp™' 72 is bounded in L 2 ((0,T) x R”) whereas in the following proposition 
we establish a bound on Vp™ in L 1 ((0,T) x R") using (16.31) . 

Proposition 6.4. There exists a constant C which does not depend on h such that 

lb/TIU 1 ([0,T];Vr 1 . 1 (n)) < CT 


for all T > 0. 

Proof. According to (16.31) . we have 


hf \V(P(p k h +1 ))\dx 
J n 


< W 2 (p k ,p k+1 ) + hC. 


Then if we sum on k from 0 to N — 1, we get 



\V(P(ph))\dxdt 






iV-l 

J2 W 2( P k h ,p k h +1 ) + TC 

k =0 
N—l 

Nj2w?(p k h ,p k+1 ) + TC 

k -0 

CT, 


because of ( 1 ( 01 ) . 

If F{x) = xlog(x) then P'(x) = 1 and if F satisfies (12.51) . then F"{x) ^ Cx m ~ 2 and P'(x) = xF"(x) ^ 
Cx m ~ 1 . In both cases, we have P'{x) ^ Cx m ~ 1 (with m = 1 for xlog(x)). So 



|V(P(pfe))| dxdt 


f [ P'(p h )\Vp h \dxdt > C f [ pn-^Vphldxdt = C [ [ |VpJ?| dxdt, 

Jo Jn Jo Jn Jo Jn 


Which proves the proposition. 


□ 


Now we introduce Q : L m (H) [0, +oo] defined by 


Sip) ■■= 


\\p m \\BV(n) if P € V ac m and p m e BV(Sl), 
-f-oo otherwise. 


Proposition 6.5. Q is lower semi-continuous on L m (Cl) and its sublevels are relatively compact in L m (fl). 

Proof. First we show that Q is lower semi-continuous on L m (ff). Let p„ be a sequence which converges strongly 
to p in L m (D) with sup„ G(pn) ^ C < + 00 . Without loss of generality, we assume that p n converges to p a.e. 
Since C < + 00 , the functions p™ are uniformly bounded in BV(fl). So we know that p™ converges weakly 
in to p. But since D is smooth and bounded, the injection of BV[fl) into L 1 ( ft) is compact. We can 

deduce that p = p m and p™ converges to p m strongly in L 1 (fl). Then by lower semi-continuity of the BV- norm 
in L 1 , we obtain 

Sip) < liminf G{p n ). 

n 

Now, we have to prove that the sublevels, A c := {p £ : G{p) ^ c}, are relatively compact in L m ( Cl). 

Since i : rj £ L 1 ( Cl) 1 —> p 1 /" 1 £ L m {VL) is continuous, we just have to prove that B c := {p = p m : p £ A c } is 
relatively compact in L 1 (D). So to conclude the proof, it is enough to notice that B c is a bounded subset of 
BV(fA) and that the injection of BV(fA) into L 1 (D) is compact. □ 
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Now we can apply Rossi-Savare theorem (theorem 15.311 to have the strong convergence in L m (]0, T[xf2) of 
Ph to p and then we find the strong convergence in L 1 (]0,T[xf2) of P{ph) to P(p), for all T > 0, using the fact 
that P is controlled by x m (12.51) and Krasnoselskii theorem (see [TO], chapter 2). 

Moreover, since 



\V{P(ph))\dxdt < CT, 


we have 

V(P(p h ))dxdt p in M n ([ 0, T] x 12), (6.7) 


i.e 


[ £-'V(P(ph))dxdt —>• [ f dp, 
J n J o in 


for all £ £ Cb([0,T] x 12) (this means that we do not require £ to vanish on <912). But since P(ph) converges 
strongly to P(p) in L 1 ([0,T] x 12), p = V(P(p)). 

To conclude, we pass to the limit in (16.41) and theorem [(Til follows. 


7 Uniqueness of solutions. 


In this section, we prove uniqueness result if 12 is a convex set. The convexity assumption is important because 
uniqueness arises from a displacement convexity argument. 

Without loss of generality, we focus here on internal energy defined on the subset of probability densities 
with finite second moment Vf c {Ll) and given by 

Vpev$°(n), Hp)~ [ F(p(x))dx, 

J n 

with F : R + —> R a convex function of class C 2 ((0, +oo)) with F( 0) = 0. We recall that for all p,p in Vf 0 , 
there exists (see for example [6]. |18] .[19]) a unique optimal transport map T between p and p such that 

W$(p,p)=[ \x-T{x)\ 2 p{x)dx. 

J f2xf2 

The McCann’s interpolation is defined by T t := Id + t(T — Id) for any t £ [0,1]. Then the curve t £ [0,1] i-» pt, 
with p t := Tt#P, is the Wasserstein geodesic between p and p ( [T%]. [lT?] . fZ \ ). 

An internal energy T is said displacement convex if 


t £ [0,1] i-t’ F(pt) is convex. 

Moreover, we say that F : [0, +oo) —>• R satisfy McCann’s condition if 

x £ (0,+oo) i-» x n F(x~ n ) is convex nonincreasing. (7-1) 

McCann showed in m that if F satisfy urn then T is displacement convex. 

Now we will state a general uniqueness argument based on geodesic convexity. This result has been already 
proved in [5] in the flat-torus case and the proof is the same in our case. 

Theorem 7.1. Assume that V) satisfy (El and (El and Fi £ satisfy ED- Let p 1 := {p\p\) and 
p 2 := (p 2 ,..., p 2 ) two weak solutions of (11.11) or (16.11) with initial conditions p\(f), •) = p\ 0 and p 2 (0, •) = p 2 0 . 
If for all T < +oo, 

nT 1 pT i 

/ J2 W v lt\\L 2 (p\)dt+ / ^2\\v 2 t \\ L 2 {p 2 )dt<+oo, (7.2) 

./0 , , Jo ;_ 1 


with, for j £ {1,2}, 


then for every t G [0, T], 


i —1 ^ u i=l 


j V? iKt) V7 TA r n 

v it : = - 1 - vv ilPt], 

Pi,t 


W 2 (p\,p 2 )^e iCt W 2 {pl,p 2 ). 

In particular, we have uniqueness for the Cauchy problems EH) and EH- 
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In the following proposition, we will prove that assumption (El holds if fZ is a smooth bounded convex 
subset of 1" or if fZ = R". 

Proposition 7.2. Let p := (p\,..., pi) be a weak solution of (11.11) obtained with the previous semi-implicit 
JKO scheme. Then pi satisfies El for all i G [1, Z]. 

Proof. We do not separate the cases where fZ is a bounded set or is R ra . We split the proof in two parts. First, 
we show that (17.21) is satisfied by p,;j, defined in (13.21) . Then by a l.s.c argument we will conclude the proof. 


• In the first step, we show that pi t h satisfies 


fj 

Jo Jn 


\^ F i(Pi,h) + Wi[p h \\ 2 pp h dxdt < C, 


(7.3) 


where C does not depend of h. 

By equality (16.31) and remark f6.31 we have 

VFl(p l *X 1 ) + V v i [p k h ] = 

where T% is the optimal transport map between 
p^f 1 and integrate on fZ, we find 


Tk(y) - y 
h 

PiX 1 and Pin- 


piX - a - e on 

Then if we take the square, multiply by 


J n IVF'^ 1 ) + Wfiptw 2 P k + 1 dx < ^W2(p*X\pth). 


Now using (12.41) . we get 

IVF/07+ 1 ) + ^Vi[p k h +1 }\ < |+ S7V i [p k h )\ + \S7Vi[pl) - Wi[p k h +1 }\ 

< I ^F'MXX) + WiiPh] I + CW 2 (p k +\p k h ) 

So we have 

j n WF'MX 1 ) + Wi[^ +1 ]| 2 ^,r dx < C (±W 2 {p k X\plh ) + W 2 (p k+ \p k )^ . 

Then using (13.911 . we finally get 


n 

l o Jn 


v*;w)+ vnp h }\ 2 pi,h dxdt = hjx [ WF'ip'fx 1 ) + wip^YpiX 1 dx 

k=0 

< cllYXwiip^Xpt^+i) 

\ k=0 / 


< c. 


• To conclude, we have to pass to the limit in (17.31) . First, we claim that \/Pi(pi,h) converges to SlPfipf) in 
Af ra ([0, T\ x fZ). In a bounded set, this has been proved in (16.7(1 . In R n thanks to the previous step, we 
have 

[ f \VPi(pi lh )\dt 
Jo Ir" 


which gives the result because Pfipi^h) strongly converges in L 1 ([0,T] x R n ) to Pfipi). 


[ f \VFl(p ith )\p i>h dxdt 
Jo J R" 

^ T ^(|VF/(p iife )| 2 + lK, 


a 
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Let 0 : R" +1 —>• Ru{+oo} defined by 


ip{r, m) 


- if (r,m) G]0,+oo[x R™, 
0 if (r, m) = (0,0), 

+oo otrherwise, 


as in [3]. And define d' : A4((0,P) x fi) x A4 n ((0,T) x SI) —>■ RU{+oo}, as in [3, by 


*(P,E) 


f* / 0 if(dp/dC, dE/dC) dxdt if p ^ 0, 
+oo otrherwise, 


where da/dC is Radon-Nikodym derivative of cr with respect to £i[o,T]xn- We can remark that since 
■0(0, m) = +oo for any m ^ 0, we have 


4' (p, E) < +oo => E <C p. 


With this definition, we can rewrite (17.31) as 

®{Pi,h,VPi(p it h)+ Wi[p h ]p it h) = [ [ | VF!(pi th ) + VVi[p h ]\ 2 p ith dxdt ^C, 

Jo Jn 

which, in particular, implies that V Pi(pi.h) <C Pi t h <C P|[o,T]xn- 

Moreover, according to [3], 4' is lower semicontinuous on A4([0,P] x SI) x _A/f n ([0,P] x SI). So, it holds 


VP, ; (/ 0 i) +Wi[p\pi) s£ liminf 4>(pi,k, VPi(p it h) + VV^P/Jp*,/,) < C, 

h \0 

which imply S7 Pfipf) <C Pi -C P|[o,T]xn an d conclude the proof because 

/ 

Jn 

= 4'(p l , yPi{pi) + Wi[p\pi) 

< c. 


VPi(pi) 


+ VK:[p] 


Pi dxdt = 


WPi(Pi) + ^Vi[p} Pi \ 


■ dxdt 



□ 
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